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between earlier results [2,3] and those presented 
here is almost entirely due to the change from 
eq. (1) to eqs. (2) and (3). 

In figure 7 we show a comparison of results 
for a s based on lattice and perturbative QCD. 
Since the analysis of ref. [5] is somewhat different 
from ours, the results labeled "NRQCD" in fig- 
ure 7 are different from the values we quote but 
agree within errors. The quarkonium results are 
in good agreement with those from other experi- 
ments, with comparable uncertainties. 



n — | — i — i — i — i — | — i — i — i — i — | — i — i — i — i — | — i — i — i — i — | — 

Re+e " 

— 9 — e + e (ev. shapes) 
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— e — LEP (ev. shapes) 
o LEP (resummed) 

6 SLD (jet rates) 

-e- Dig (u) 
-e- Dig (fx) 
-e- CLEO (T decays) 

— b— bb 
-H- NRQCD 

B CC 

□ NRQCD 

j I i i i i I i i i i I i i i i I i i i i I 

0.10 0.12 0.14 0.16 0.18 

« (5 W m z) 

Figure 7. a |jg-( m z) from lattice QCD in compar- 
ison with perturbative results. For recent refer- 
ences see ref. [1]. 



5. CONCLUSIONS AND OUTLOOK 

An important goal of lattice QCD is the de- 
termination of Standard Model parameters from 
experimental input. Among the first results that 



include all systematic errors are determinations 
of a s from quarkonium spectroscopy. This is an 
important achievement, even though the calcula- 
tions are not entirely based on first principles. 

The phenomenological correction that accounts 
for the systematic effects from the omitted sea 
quarks has been tested in two ways. The consis- 
tency of the results from bb and cc spectra is a 
phenomenological test of the correction. A first 
principles test with dynamical fermions was pro- 
vided by the KEK collaboration [10], albeit with 
large errors. 

As already mentioned in the introduction, phe- 
nomenological corrections are a necessary evil 
that enter most coupling constant determina- 
tions. However, in the long run, lattice QCD will 
give us a truly first principles determination of a s 
with an accuracy that cannot be matched by any 
other method. At that point, the MS standard 
should be replaced with a new one, a coupling 
that can be computed both in lattice QCD and 
perturbative QCD. 
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Figure 5. The and 3,4 flavor running of ay from 
a matching scale at q = 700 MeV to q = 5 GeV. 

Lepage e< a/. [5] suggest to estimate q* as the ex- 
pectation value of the potential. One could also 
use the scale at which the slopes of the nj = and 
3 flavor potentials match. These different pro- 
cedures lead to estimates for q* consistent with 
350 - 650 MeV for cc states and 600 - 1000 for bb 
states. Note that the typical momenta in the p- 
wave states are more infrared than the momenta 
in s-wave states. The uncertainty in q* leads to 
an estimate of the uncertainty of the sea quark 
correction. 

Before we go on to calculate the correction, a 
few words on the two-loop evolution of the cou- 
pling are in order. Lepage et al. [5] show that ay 
runs with the 2-loop /3-function down to momenta 
~ 700 GeV (see figure 1 in ref. [5]). Something 
similar was also observed in ref. [12]. We thus 
conclude that 2-loop evolution is probably okay 
for bb. We find a 20% correction, with a 5% un- 
certainty, to ttf^g - at 5 GeV. However, close to the 
Landau pole, the perturbative coupling diverges 
and nonperturbative effects become important. 
Consequently, at the low end of the cc matching 
scales the perturbative correction becomes unsta- 
ble and has to regulated with a potential model. 
We find a 26% correction with a 8% uncertainty. 

A first principles test of the correction has now 
been performed for the first time by the KEK 
group, as shown in figure 2 of ref. [10]. The per- 
turbative evolution of and 2 flavor couplings to 



Table 3 

Summary of uncertainties in the a s determination 
at 5 GeV. 



uncertainty 


bb (%) 


cc (%) 


stat., a 2 


1 - 3 


1 - 3 


perturbative 


5 


5 


n f = ~> n f 


= 3,4 5 


8 


total 


8 


10 


total at mz 


5 


8 



the charmonium scale gives consistent results, al- 
beit with large errors. 

(5) 

Figure 6 shows the final result for a-^r, evolved 
to the Z mass scale, for the different cases (with 
statistical and systematic errors included). With 
the exception of the KEK result, the determina- 
tion is totally dominated by the systematic errors 
from the quenched approximation and perturba- 
tion theory, leaving the other lattice errors ne- 
glegibly small in comparison. Thus, the quenched 
results are in much better agreement than the er- 
rorbars shown. 

0.14 r □ 

"n 0.12 — T — 

6 ; i i i T = 

X o.io - 1 -J- - 

Is : : 

2 0.08 ~— —_ 

S : bb cc 

0.06 t =1 

Figure 6. a ^^( m z) m comparison; Fermilab (□: 
bb, O cc), NRQCD (V: bb, x: cc), KEK (*), 
UK(NR)QCD (o), and UKQCD (+). 



Table 3 lists all the uncertainties that enter the 
lattice determination of a s . We quote as our final 
result the value obtained from the bb spectrum: 

a^I(m z ) = 0.108 ± 0.006 . (6) 

The cc result, a ^^( m z) = 0.110(8), is in good 
agreement with our central result. The difference 
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Figure 3. The nonperturbative coupling of ref . [5] in comparison with perturbative predictions. 



cc that was observed in a -1 is still present. Also, 
the 2 flavor coupling [10] is about 10% bigger than 
the quenched results. 



0.18 
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Figure 4. a ^-(5 GeV) vs. a 2 . Fermilab (□: 
bb, O: cc), NRQCD (V: bb, x: cc), KEK (*), 
UK(NR)QCD (o), and UKQCD (+). 



In the next section we will show how, after 
proper inclusion of sea quark effects, everything 
falls into place. 



4. SEA QUARK EFFECTS 

This final correction can, at this point, not be 
included entirely from first principles. However, 
it can be estimated using potential model phe- 
nomenology and perturbation theory [3,14]. 

The physical quantity that is used to determine 
a -1 sets the "matching scale" , where by construc- 
tion the coupling of the quenched lattice poten- 
tial matches the coupling of the effective 3 flavor 
potential. The difference in the perturbative (2- 
loop) evolution of the and 3,4 flavor couplings 
then causes the quenched coupling at short dis- 
tances to be too small. This is illustrated in fig- 
ure 5. 

Different "matching scales" for bb and cc then 
lead to differences in the quenched coupling at 
short distances, as observed earlier. The match- 
ing scale is, in effect, the average gluon momen- 
tum transfer in the 'onium state in question. It 
can be estimated as the expectation value of the 
gluon momentum using potential model or lattice 
wave functions: 



{q)is,iP,2S,. 



(5) 



3 




5.6 



5.8 



6.0 



6.2 



Figure 2. Summary of a 1 vs. (3 = 6/<7q from 
the following groups: Fermilab [1] (□: 66, O: 
cc), NRQCD [2] (V: 66, x: cc), KEK [3] (*), 
UK(NR)QCD [4] (o), and UKQCD [5] (+). 



3. RENORMALIZED COUPLING 

Within the framework of Lattice QCD the con- 
version from the bare to a renormalized coupling 
can, in principle, be made nonperturbatively. 
One can, for example, define a renormalized 
coupling from the nonperturbatively computed 
heavy quark potential (see the talk by Schilling 
[11]). In [12] a renormalized coupling is defined 
nonperturbatively through the Schrodinger func- 
tional; its running is computed using a finite size 
scaling technique. Of course, as always, the asso- 
ciated lattice spacing errors have to be carefully 
removed. 

In ref. [13] it was shown that perturbation the- 
ory describes all short distance lattice quantities 
considered there up to a few percent provided a 
renormalized coupling is used. In ref. [2,3] the 
MS coupling was determined from the bare lattice 
coupling using the mean field formula proposed in 
ref. [13]: 



An alternative is to define a renormalized cou- 
pling through short distance lattice quantities, 
like small Wilson loops or Creutz ratios. For ex- 
ample, the coupling defined from the the plaque- 
tte, ap = — 3 In (TWp)/47r, can be expanded in 
ay by [13]: 



a P = a v (q)[l - l.l9a v (q) + 0{a v ) 



or, equivalently 
a v 1 (q) = a] 



l.l9 + 0(a v ) 



(2) 



(3) 



where q = 3.41/a is chosen according to the pre- 
scription in ref. [13]. The MS and V couplings 
are related by: 



\q) + 0.822 



(4) 



Air{TrU P )/gl + 0.309 



(1) 



The size of higher order corrections associated 
with the above defined coupling constants can be 
tested by comparing perturbative predictions for 
short distance lattice quantites with nonperturba- 
tive results. It was shown in ref. [13] that, while 
the mean field coupling (1) accounts for the bulk 
of the corrections, ay from the plaquette yields 
more accurate results. We will use both eqs. (2) 
and (3) for the analysis presented here, and aver- 
age over the results. 

In figure 3, the nonperturbative coupling from 
ref. [12] is compared to perturbative predictions 
for this coupling using eq. (3) and the 1-loop re- 
lations given in ref. [12]. Similarly, in ref. [5] the 
reliability of ay using eq. (2) is estimated, from 
the spread observed in small Wilson loops and 
Creutz ratios up to size 3 compared to the plaque- 
tte. The next to next to leading order corrections 
for these quantities have been calculated numer- 
ically. We estimate a 5% uncertainty associated 
with perturbation theory. 

Figure 4 shows the zero-flavor MS coupling (us- 
ing eq. (2)) as a function of a 2 for 66 and cc for the 
different groups in comparison (with statistical 
errorbars only). The Fermilab results alone indi- 
cate that residual lattice spacing artifacts (higher 
than 0(a)) are too small to resolve with the sta- 
tistical errors present. The scatter between re- 
sults from the different groups, using either 66 or 
cc states, is within the range of statistical errors. 
The small difference in a s from 66 as compared to 
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Table 1 



Comparison of various determinations of a s . 



process 


observable 


theory 


caveats 




had. event shapes, . . . 


pert. QCD (nlo, resummed) 


hadronization 


e + e~ 


R e + e- , Rz 


pert. QCD (nnlo) 






R T 


pert. QCD (nnlo), analyticity 


OPE, QCD sum rules 




T decays 


pert. QCD (nlo) 


relativistic corrections 


DIS 


^3 


pert. QCD (nlo) 


higher twist 




F 2 


pert. QCD (nlo) 


higher twist, g(x,Q 2 ) 


pp 


pp —> W+ jets 


pert. QCD (nlo) 


large pert, correction 


bb, cc 


spin averaged splittings 


lattice QCD 


sea quark effects 



2. SETTING THE SCALE 

In setting the scale there are two sources of sys- 
tematic error to be considered, apart, of course, 
from sea quark effects. 

The leading lattice spacing errors of the rel- 
ativistic Wilson action are spin dependent and 
therefore do not contribute to spin-averaged 
quantities. Likewise, for the non-relativistic 
action, spin-dependent interactions are added 
as higher dimensional operators; spin-averaged 
quantites do not depend on the additional cou- 
plings. 

The scale, a -1 , depends in principle also on 
the other parameter which enters the calculation, 
the quark mass [4]. Table 2 shows that the ex- 
perimentally observed mass dependence for spin- 
averaged splittings between the cc and bb systems 
is mild. 

Figure 1 shows the mass dependence of the 
IP-IS splitting in a quenched lattice calculation 
by the Fermilab group. Note that the splitting 
is about 10% smaller for bb than for cc. Also 
shown in figure 1 is the (in-)dependence of the 
IP-IS splitting on the improvement term. It is 
thus clear that spin-averaged splittings give the 
most accurate determination of the scale, a -1 . 
The NRQCD collaboration fits to the entire spin- 
averaged spectrum to determine a -1 [4,5]. 

A summary of all the results for a -1 is given 
in figure 2 as a function of the gauge coupling 
(j3 = 6/<7q)- The comparison of a -1 from the bb 
and cc spectra shows a small scale dependence 
which is already apparent in figure 1. This am- 
biguity can only be resolved with the proper in- 
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Figure 1. The IP-IS splitting vs. kinetic mass on 
the 24 3 x 48 lattice at j3 = 6.1. x : Wilson (c = 0) 
action; □: Improved (c = 1.4) action. 



Table 2 

Spin-averaged splittings in the T and J/ip sys- 
tems in comparison. 



quantity 


cc (MeV) 


bb (MeV) 


m(lP - IS) 


456.8 


452 


m(2S - IS) 


596 


563 


m(2P - IP) 




359.7 



elusion of sea quark effects. In order to see how 
consistent the different results for a -1 are with 
each other, we will apply in the following the same 
analysis in converting from a -1 and a\ at = 47r 
to «gl(5 GeV). 
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The status of determinations of a s from quarkonia using lattice QCD is reviewed. We compare the results with 
those obtained from perturbative QCD. 



1. INTRODUCTION 

At present, the QCD coupling, a s , is deter- 
mined from many different experiments, per- 
formed at energies ranging from a few to 90 
GeV [1]. In most cases perturbation theory is 
used to extract a s from the experimental infor- 
mation. Experimental and theoretical progress 
over the last few years has made these deter- 
minations increasingly precise. However, all de- 
terminations, including those based on lattice 
QCD, rely on phenomenologically estimated cor- 
rections and uncertainties, from nonperturbative 
effects. Table 1 (adapted from S. Bethke [1]) lists 
the caveats associated with each determination. 
These effects will eventually (or already do) limit 
the accuracy of the coupling constant determi- 
nation. In the case of the lattice determination 
the limiting uncertainty comes from the (total or 
partial) omission of sea quarks in numerical sim- 
ulations. We review the progress that has been 
made since the first determinations of a s from the 
lattice [2,3]. 

The lattice determination of the strong cou- 
pling, a s , requires three ingredients: 

1. The experimental input to the lattice de- 
termination is usually a mass or mass splitting, 
from which by comparison with the correspond- 
ing lattice quantity the scale, a -1 , is determined 
in physical units. For this purpose, one should 
identify quantities that are insensitive to lattice 
errors. 

2. In the definition of a renormalized coupling, 



systematic uncertainties should be controllable, 
and at short distances, its (perturbative) relation 
to other conventional definitions calculable. 

3. Calculations that properly include all sea 
quark effects do not yet exist. If we want to 
make contact with the "real world" , these effects 
have to be estimated phenomenologically. One 
should therefore in "setting the scale" consider 
physical systems where the reliability of the phe- 
nomenological correction can be quantitatively 
estimated. 

As has been pointed out some time ago by Lep- 
age [2], quarkonia are the ideal systems for this 
program, since systematic errors can be analyzed 
using potential models. 

Several groups presented results for quarko- 
nium spectroscopy, which are reviewed by Davies 
[4]. The NRQCD collaboration [5,6] has re- 
sults for the bb and cc spectra using a nonrel- 
ativistic QCD action at next to leading order. 
The leading order NRQCD action is used by the 
UK(NR)QCD collaboration [7] who presented re- 
sults for the spin-averaged bb spectrum. The Fer- 
milab group [3,8] has results for both the bb and cc 
spectra, using an O(a) improved Wilson action. 
The UKQCD collaboration [9], using a similar ac- 
tion, looked at the cc spectrum. The first calcu- 
lation of the cc spectrum with dynamical (nj = 2 
staggered) fermions was presented by Onogi [10] 
for the KEK collaboration; they looked at the 1P- 
1S splitting using Wilson and staggered valence 
fermions. 
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